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Abstract Phylogenetic relationships among the populations of Erebia ligea on the Japanese archipelago were inferred 
from mitochondrial DNA sequences. The results show a clear division into three lineages on the Sakhalin, Hokkaido and 
Honshu islands, and also suggest that these lineages have been reproductively isolated from each other since before the 
last ice age. Further fragmentations were found within these lineages: two genetic lineages on Hokkaido (one found at 
higher altitudes in the Daisetsuzan Mts. and the other at lower altitudes around some mountain bases) and two on Honshu 
(one found in the Hida, Yatsugatake and Hakusan Mts. and the other in the the Akaishi and Kiso Mts.). Nested clade 
phylogeographical analysis (NCPA) points to a history of E. ligea in the Japanese archipelago which was punctuated by 
several repeats of allopatric fragmentation and restricted gene flow. The existence of a secondary contact zone shows the 
same scenario. In the Hida Mts., different haplotypes were detected in the northern and the southern regions, with a 
secondary contact zone found in the vicinity of the central region. The ancient population was fragmented, reproductively 
isolated and then differentiated into several populations with distinct haplotypes. Distributions of the two populations 
then expanded allowing a secondary contact of the two populations, and sporadic distributions were subsequently formed. 
Repetitive cycles of expansion and reduction in distribution within the Japanese archipelago in response to the glacial 
cycles of the Quaternary period are suggested. 
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Introduction (Fig. 1). The Tartar Strait and Soya Strait are shallow, and 
they were formed less than 11,000 years ago (Ohshima, 
2000). On the other hand, the Tsugaru Strait and Tsushima 
Strait are deep. Honshu separated from the Korean peninsula 
and Hokkaido, and became an island during the interglacial 
age after the Riss glaciation about 150,000 years ago. It 
is believed that Honshu has not been connected to the 
continent since then (Ohshima, 2000). Sakhalin and 
Hokkaido had been connected to the continent almost 
continuously during a period of 65,000 years from about 
1.1 Phylogeographic history of alpine organisms on 75,000 years ago to about 10,000 years ago. They had 
the Japanese archipelago been parts of the continent for longer periods of time than 
the periods when they were isolated islands. Thus, during 
the glacial period geographic conditions allowed organisms 
to migrate southward from the continent through Sakhalin 
to Hokkaido (Ono, 1990). On the other hand, alpine 
organisms that now have an isolated distribution in the 
alpine region of Honshu had expanded their distribution 
prior to about 150,000 years ago from the continent, and 
they are the species that were affected by climate changes 
since then during the interglacial period and the last glacial 
period. 


Alpine organisms are most notably affected by global 
warming, and therefore reduction of warming by human 
acts and conservation of these organisms are necessary. 
In this paper, we focus on alpine butterflies as a representative 
and assess the history of distribution of populations in the 
Japanese archipelago during the Quaternary period and 
their genetic relationships with the continental populations. 


There were periods when the Japanese archipelago was 
connected to the Asian continent during glacial ages and 
served as a refugium for the animals and plants of 
continental origin. During glacial ages, many of the arctic 
and subarctic animals and plants migrated into the Japanese 
archipelago. During the following interglacial age, ranges 
of those animals and plants diminished and many of them 
disappeared. There have been many paleontological studies 
of fauna based on fossils (Kawamura, 1998), and of flora 
based on pollen and wooden fossils (Suzuki and Takeuti, 
1989). High altitude areas of central Honshu and Hokkaido in 
the Japanese archipelago have had at least two periods 
when glaciers developed (Iwata, 2003). Reconstruction 
studies of ancient flora by pollen analysis (for example, 


The major islands that constitute the Japanese archipelago 
are Sakhalin, Hokkaido, Honshu, Shikoku, and Kyushu 
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Fig. 1. The Japanese Archipelago. 
Major islands and straits are indicated. The small 
rectangle shows the location of Hida Mts. 





Fig. 2. Erebia ligea takanonis (photo courtesy of M. Fukumoto) 


Hoshino, 1998) suggest that low altitude areas on Hokkaido 
were the southern limit of taiga during the last glacial 
period, and that there were forests richer with Larix gmelinii 
than those in the present northern Sakhalin. In the eastern 
and northern parts of Hokkaido, the distribution of taiga 
with Larix gmelinii as its main plant were divided in such 
a way that coniferous forests, grasslands, wet vegetation 
and alpine tundra formed segregated niches according to 
altitude and/or geographical features. Thus the areas 
functioned as refugia for many organisms in the Far East 
(Ono, 1990). 


The limitations on regional biogeographic approaches or 
studies based on-the-degree-of-specialization-of external- 
characters for understanding immigration routes of alpine 
butterflies into the Japanese archipelago, and varying 
patterns of their distribution on the islands, have been 
pointed out (Hiura, 1977). Recent studies based on a new 
viewpoint provided by a molecular phylogenetic approach, 
circumvent these limitations. Some examples are as follows: 
studies on the alpine butterflies Erebia neriene ssp. 
(formerly treated as Erebia niphonica) (Nakatani et al., 
2007a), Oeneis norna ssp. (Usami et al., 2011), and also 
on some alpine plants (Fujii and Senni, 2006). 


The process of expansion and contraction of the ranges 
of organisms due to climate changes in the Quaternary is 
well known as a process of north-south migration on the 
Japanese archipelago. According to our studies on alpine 
butterflies, small populations of both Erebia neriene ssp. 
and Oeneis norna ssp. survived during the interglacial 
period within more than one refugium on the Japanese 
archipelago. As a result of isolation in each refugium 
during the warm periods, genetically differentiated 
populations have survived till now, forming multiple 
lineages by expanding and contracting their distribution 
areas as they adapted to later glacial cycles (the multiple 
refugia model). On the other hand, some species of alpine 
plants have two lineages in Honshu, northern and southern, 
and the boundary is in the northeast of the island. As the 
primary cause of such a distribution character, the following 
scenario has been suggested. The original population 
remained in the mountainous areas of central Honshu as 
a refugium during the interglacial period and later became 
the southern lineage. A new lineage then migrated into 
the Japanese archipelago during the next glacial period 
and formed the northern lineage (the single refugium 
model). 


1.2 Systematics and ecology of the subject species 


Erebia ligea (Fig. 2) is widely distributed in the mid to 
high latitude regions of the Eurasian continent, from the 
European Alps and Scandinavia in the west through the 
Ural mountains to the Kamchatka Peninsula in the east. 
Populations are found in the Japanese archipelago on 
Sakhalin, Hokkaido and Honshu islands, and are recognized 
as separate subspecies (ssp. rishirizana from Sakhalin and 
Hokkaido, ssp. takanonis from Honshu). In Honshu, this 
species inhabits grasslands in coniferous forests where 
species such as Poa and/or Carex grow, or grasslands just 
above the timberline, and the distribution ranges are found 
from 1,700-2,500 m alt. Although E. ligea frequently 
coexists with E. neriene in Honshu, E. ligea has more 
limited distribution areas. The upper elevation limit 
particularly is lower, while the lower elevation limit is 
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higher, leaving only a narrower zone of suitable elevation. 
In Hokkaido, E. ligea is found at a higher elevation in the 
Daisetsuzan Mts., and also at the foot of some mountains 
approximately 1,000 m in height, often along woodland 
paths. An isolated population is also found at the summit 
of Rishiri Island, northern Hokkaido. This kind of 
distribution pattern, present, though sporadically, in wide 
areas is desirable for a study of phylogeographic history 
(DeChaine and Martin, 2004; Nakatani et al., 2007a). 


2. Materials and methods 
2.1 Specimens and collection of molecular information 


The authors collected 84 individuals from the Japanese 
archipelago, and also 12 individuals in Europe, Mongolia, 
and Russia for reference. The bodies of the sample 
specimens without wings were preserved in 99% ethanol 
and kept at 4°C until DNA was extracted. 


For genetic markers, we used the NADH dehydrogenase 
subunit 5 (NDS) and cytochrome oxidase subunit 1 (CO7) 
genes in mitochondrial DNA. DNA was extracted from 
the thoracic flight muscles or legs by DNeasy Tissue 
Extraction Kit (Qiagen, Hilden, Germany). PCR was 
conducted in a 10ul reaction system by adding DNA 
(approximately 40ng) with reaction constitution as 
recommended by Applied Biosystems. Thermal cycling 
profile was one step at 95°C for 10 min; 25 cycles at 95°C 
for 30 sec, 45°C for 30 sec and 72°C for 1 min 30 sec; one 
extension step at 72°C for 8 min 30 sec. After amplification, 
the PCR products were purified with ExoSap-IT (USB, 
Cleveland, OH, USA). The nucleotide sequences of the 
amplified DNA fragments were determined by the direct 
sequencing method with BigDye Terminator v1.1 Cycle 
Sequencing Kit (ABI, Weiterstadt, Germany) on an ABI 
Genetic Analyzer. As primers for PCR and direct 
sequencing, we used V1 (5'-CCT GTT TCT GCT TTA 
GTT CA-3), C2 (5'-ATC YTT WGA ATA AAA YCC 
AGC-3') (Yagi et al., 1999) for ND5, and mtD6 (C1-J- 
1718; 5'- GGA GGA TTT GGA AAT TGA TTA GTT 
CC-3'), Nancy (C1-N-2191; 5'-CCC GGT AAA ATT AAA 
ATA TAA ACT TC-3') (Simon et al., 1994) for CO1. We 
conducted fragment comparison of 432 bp for NDS, 465 
bp for CO1, 897 bp in total. 


2.2 Sequence alignments and character statistics 


Mitochondrial ND5 and CO/ sequences were edited and 
aligned using SeqScape V. 2.11 (ABI, Weiterstadt, 
Germany). Base frequency homogeneity was tested 
separately for each dataset with the y?-test using PAUP* 
4.0b10 (Swofford, 2002). The X 人 test did not reject the 
hypothesis of homogeneity of nucleotide frequencies in 
every pair of taxa (P = 1.00). To test conflicts in phylogenetic 


signal among each dataset we conducted an incongruence 
length differences (ILD) test (Farris et al., 1994) in PAUP* 
4.0b10 with heuristic searches with tree bisection and 
reconnection (TBR) and 100 random addition replicates 
for each. The ILD test revealed no conflict between ND5 
and CO1 (P=0.58). 


2.3 Nested clade phylogeographical analysis 


Nested clade phylogeographical analysis (NCPA) 
(Templeton et al., 1995; Templeton, 1998, 2004, 2008) 
can infer past events that affected the population constituents 
of a species such as past fragmentation, restricted gene 
flow and range expansion by combining haplotype 
phylogeny and geographical information on its distribution. 
We inferred topology among haplotypes by using statistical 
parsimony methods (Templeton et al., 1992) integrated 
into the TCS v. 1.21 (Clement et al., 2000). The network 
is connected with 95% reliability that multiple transfers 
do not occur at any locus. The relationships between the 
inferred phylogeny and geographic position can be 
evaluated by comparison of two distances of the within 
clade distance (Dc) and the nested clade distance (Dn), 
and by the interior-tip contrasts (I-T). The significance of 
the distances defined by NCPA is evaluated using a=0.05. 
These evaluations were calculated by the GeoDis program 
(Posada et al., 2000). When the evaluation is significantly 
large (L) or significantly small (S), they indicate that 
haplotype/clade distribution is significantly widespread 
or significantly restricted. Significant I-T contrast for both 
Dc and Dn values denotes a displaced (L) or restricted (S) 
distribution of lower versus higher haplotypes/clades. 


3. Results 


3.1 Genetic difference among the continent, Sakhalin, 
Hokkaido, and Honshu populations 


Distribution frequency of haplotypes and sample sites are 
shown in Table 1. As a result, 17 haplotypes from the 
Japanese archipelago (including Sakhalin) and 7 from the 
continent were detected. Topology of the relationships 
among haplotypes was inferred using a statistical parsimony 
method (Fig. 3). Due to the large number of mutational 
steps separating the continental and the Japanese archipelago 
haplotypes, they could not be joined in the network with 
greater than 95% confidence. C り / (465 bp) fragment 
comparison indicated that there is about 2.9% (uncorrected 
pair wise distance) differences of substitutions between 
major haplotypes of the continental and the Japanese 
archipelago populations. Also, between the populations 
in Honshu and Hokkaido, or those in Hokkaido and 
Sakhalin, there is about 0.8% and 1.6% difference 
respectively. Within the continent, the difference between 





NII-Electronic Library Service 


The Lepidopterological Society of Japan 





Phylogeography of Erebia ligea in the Japanese archipelago 207 


Table 1. Distribution of the haplotypes (NDS plus CO/) of Erebia ligea from each sampling range. 
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populations taken in Europe and those in Mongolia is only 
0.2%. 


3.2 Nested clade phylogeographical analysis (NCPA) 
on Honshu populations 


It is clear that the populations of Sakhalin, Hokkaido, and 
Honshu were separated from each other by past allopatric 
fragmentation. So we studied the distributional transition 
in Honshu. All haplotypes from Honshu were grouped 
into a single 3rd level clade (Fig. 5). NCPA showed a 
significant association between genetic structure and 
geographic distribution (Fig. 6). There are 5 clades where 
a significant tendency (P<0.05) of genealogy and geographic 
location was suggested by NCPA (Table 2). 


In the application of NCPA, the age of a clade is the 
maximum of the ages of the smaller clades nested within 
it (Templeton et al., 1995). Hence, the older clades are in 


general the larger clades with more nested series of clades. 
Three clades, 1-1, 2-1, and 3-1, were within a nested series 
of clades, and also appeared to be statistically significant. 


(1) Early fragmentation 


After migrating into the Japanese archipelago, and spreading 
out into areas that include the present distribution ranges, 
they dispersed into three disjunct refugia during warm 
interglacials, and divided into three lineages, the Sakhalin, 
Hokkaido, and Honshu lineages (allopatric fragmentation). 


(2) The second fragmentation 


A panmictic population at an early glacial period divided 
into three separate populations of the northern-1 lineage 
(southern Hida Mts. and Yatsugatake Mts.; clade 2-1), 
northern-2 lineage (northern Hida Mts. and Hakusan Mts.; 
clade 2-2), and the southern lineage (Akaishi Mts. and 
Kiso Mts.; clade 2-3) as a result of isolation in three refugia 
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Fig. 3. Statistical parsimony network. 
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Circles denote haplotypes, the size of the circle reflects the wideness of distribution ranges in 2 
levels, small solid circles denote missing or theoretical haplotypes, and a node represents a single 
mutational step. The pie chart shows the proportions of regions with the same haplotype. Uncolored 


haplotype was detected in one region only. 


during the following warm interglacial period (past 
fragmentation, clade 3-1). 


(3) Isolation by distance 


The northern-1 lineage (clade 2-1), which is a 2nd level 
clade belonging to clade 3-1, is inferred to have two 
separate populations, the northern Hida Mts. population 
(clade 1-2), and southern Hida Mts. and Yatsugatake Mts. 
population (clade 1-1), by restricted gene flow, particularly 
gene flow restricted by isolation-by-distance. 


(4) Fragmentation and habitat islanding 


The southern Hida Mts. and Yatsugatake Mts. population 
(clade 1-1) is fragmented into the Hida Mts. and Yatsugatake 
Mts. populations. 


4. Discussion 
4.1 Distributional transitions in the Japanese Archipelago 


Based on the inference key of NCPA (Dr Posada’s web 
site, var. Jan 06, 2011), the history of E. ligea in the 


Japanese archipelago was punctuated by several cycles of 
allopatric fragmentation and restricted gene flow. Inferences 
of historical processes made by NCPA should be viewed 
with caution, because the inference keys make qualitative 
assessments (Knowles and Maddison, 2002; Panchal and 
Beaumont, 2007; Templeton, 2008). But the presence of 
traces of secondary contacts of two populations with 
different haplotypes, found in the area close to Mt. 
Harinokidake (hereinafter referred to as “Harinoki Gap”, 
Fig. 4), northern Hida Mts., is direct evidence that the 
Erebia ligea populations underwent repeated fragmentation 
and dispersal as a result of glacial cycles. It suggests that 
distribution ranges fragmented into separate areas like 
islands in the past warm period, in each of which a 
population with characteristic haplotypes evolved, and 
that later in the glacial period these populations expanded 
and came into secondary contact with each other, following 
which, in the next warm period, the distribution ranges 
have again become like islands. These traces of secondary 
contact of different haplotypes strongly support the notion 
that cycles of fragmentation and dispersal of distribution 
ranges in the Pleistocene took place in the Japanese 
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Fig. 4. Sampling sites and geographic distribution of haplotypes 
from Hida Mts. Red: HI000, blue: HI100. 


archipelago. The history is similar to that of E. neriene 
(Nakatani et al., 2007a). The existence of such traces of 
secondary contacts of different haplotypes observed for 
two distinct species in Japan strongly supports the notion 
that cycles of fragmentation and dispersal of distribution 
ranges the in Pleistocene took place in the Japanese 
archipelago. 


4.2 Population genetic structures as compared with 
those of Erebia neriene 


Conceptual distribution maps of major haplotypes and 
their phylogenetic relationships are shown for E. ligea 
(Fig. 7A) and for E. neriene (Fig. 7B). In both species, 





clear genetic differentiation can be identified between the 
populations in the Hida Mts. and Akaishi Mts. The ancestral 
population, which initially had a single genetic structure, 
was separated into two populations, one distributed in the 
Hida Mts. and another in the Akaishi Mts. through the 
process of a reduction in distribution areas during the 
subsequent warmer period. 


The genetic distances between the two lineages, originating 
from a refugium in the Hida or Akaishi Mts., are 
differentiated between the species E. ligea and neriene, 
suggesting a possibility that the two species fragmented 
around different warmer periods. In both species, the 
Akaishi lineage is closest to the northernmost population 
in the Hida lineage, allowing for the conclusion that the 
expansion routes might be similar. 


For the populations in the Hida Mts., the major difference 
between the two species is that E. neriene has a unique 
haplotype, while E. ligea is divided into two sub-lineages, 
north and south, separated with a range of sympatry in 
the Harinoki Gap. In the Akaishi Mts. conversely, E. ligea 
has a unique haplotype, and E. neriene is divided into two 
sub-lineages with a range of sympatry. Despite the Hakusan 
Mts. being geographically far distant from the Hida Mts., 
both species share the same haplotype in the two mountain 
systems (for E. ligea, the northern sub-lineage in the Hida 
Mts.), suggesting that the populations in the two mountain 
systems were, up until recently, conducting genetic 
exchanges in both species. 


E. neriene in the Yatsugatake Mts. shares the same haplotype 
as that in the Hida Mts., while the haplotype of E. ligea 
in the Yatsugatake Mts. differs from that in the southern 
sub-lineage in the Hida Mts. by only a single base 
substitution, suggesting a close relationship between the 
populations of the two regions. In contrast E. neriene in 
the Kiso Mts. shares the same haplotype as that of the 
Hida Mts., while E. ligea in the Kiso Mts. is closely related 
to the Akaishi lineage. This suggests that different scenarios 
of distribution transitions should be inferred for the 
populations of the two species in the Kiso Mts. Concerning 
the Sakhalin and Hokkaido populations, E. neriene 
populations on the two islands share the same haplotype, 
indicating recent genetic exchanges. For E. ligea, however, 
there is a large genetic difference between the populations 
on the two islands, suggesting reproductive isolation ever 
since the era before the last ice age. In contrast to the 
widespread distribution of E. neriene in the lowland of 
Hokkaido (particularly at the northern end, where this 
species is present even on terraces near the seashores), the 
distribution of E. ligea is limited to higher altitudes over 
1,000 m in central Hokkaido. The difference in population 
genetics between the two species would be supported by 
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Table 2. Inferences of historical processes 


Clade Sreps in inference key Inferred event 





ye P-value 





1-1 1-2-3-5-15-No past fragmentation 21.00 0.0030 
1-3 1-2-3-4-No restricted gene flow 46.13 0.0200 
2-1 1-2-3-4-No restricted gene How 14.60 0.0130 
2-3 1-2-3-5-15-No past fragmentation 12.00 0.0020 
3-1 1-2-3-5-16-No allopatric fragmentation 50.59 0.0000 





x°-test of geographical association of clades and inferences of biological causes for 


such association. 


The inferences were obtained following the most updated key (Jan 06, 2011). 











Fig. 5. Haplotype network and associated nested design. 
All haplotypes are nested in the 4th level clade. 


the difference in their ecological behavior described above. 


As outlined above, E. ligea and neriene seem to have 
undergone different histories of distribution transitions in 
some localities in the Japanese archipelago, although they 
are closely related species of the same genus and prefer 
similar habitats. Similar examples of closely related species 
having a different era or location in which reproductive 
isolation took place, are also observed in Europe. Erebia 
medusa is found in grasslands in coniferous forests and 
in alpine meadows above the timberline in Bulgaria. 


Schmitt et al. (2007) studied the populations from Romania 
and Bulgaria and pointed out that the species is separated 
into two distinct lineages, north and south, separated by 
the Danube valley. Despite the species having expanded 
its range after the last ice age, no genetic exchange occurred 
across the river valley. Erebia euryale is also found in 
grassland in the coniferous forest. The Romanian (the 
Carpathian mountains) and the Bulgarian (the Rila 
mountains) populations are, however, genetically very 
similar, suggesting that genetic exchanges across the 
Danube valley took place even after the last glacial period 





NII-Electronic Library Service 


The Lepidopterological Society of Japan 





Phylogeography of Erebia ligea in the Japanese archipelago 211 


Ist Step Clade 2nd Step Clade Total Cladogram 
No Dc Dn No Dc Dn No Dc Dn 
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Fig. 6. Analytical results for NCPA with significant geographic structures. 
Nesting hierarchy proceeds from the most recent Ist level clades on the left to the older 
clades on the right. Within each clade (boxed), columns indicate haplotypes, one 
hierarchical lower level clade or interior-tip contrasts, within clades distance (Dc), and 
nested clades distance (Dn). Individual haplotypes/clades with asterisks are tips. 
Significantly large and small statistical evaluations are denoted by L and S, respectively. 


(Schmitt and Haubrich, 2008). E. ligea is also found in 
grassland in coniferous forest and is often sympatrically 
distributed with E. euryale, but a strong differentiation 
between the northern and southern mountain ranges in 
Bulgaria has been shown (Nakatani et al., 2012). 


As mentioned above, for different species even sharing a 
similar habitat in a particular region in the continent, 
different genetic population structures were found, which 
had probably formed in response to the various glacial 
cycles during the Quaternary period. In the Japanese 
archipelago, different genetic structures among different 
species may also be observed derived from species-specific 
biological abilities to adapt to climatic changes in the 
Quaternary period. 


4.3 The ecological significance of the Harinoki Gap 


As mentioned previously, a sympatric distribution area of 
the two sub-lineages of E. ligea was found in the Hida 
Mts., located at the Harinoki Gap (Fig. 4). Oeneis norna 
is a true alpine butterfly found only in the alpine zone 
above 2,700 m alt., the screes and dwarf pine (Pinus 
pumila) areas. In the vicinity of the Harinoki Gap, O. 
norna is divided into a northern sub-lineage distributed 
from Mt. Yukikuradake to Mt. Nunobikiyama and a southern 
sub-lineage distributed further south of Mt. Eboshidake, 
and is absent between Mt. Jiigatake and Mt. Eboshidake 


(Usami et al., 2011). In addition, a unique lineage of E. 
neriene found in the Hida Mts. is spreading over the 
Harinoki Gap. 


The distribution patterns of these three alpine butterflies 
around the Harinoki Gap can well be understood mainly 
as the reflection of species-specific fitness for certain 
altitudes (Fig. 8). Observing the geography and vegetation 
of the Harinoki Gap (Fig. 9), a large portion of the ridge 
is covered with coniferous and broad-leaved forests (Betula 
ermanii) due to the low altitudes of 2,500-2,600 m; the 
western side of the ridge is windward to the winter monsoon 
causing the forest to grow up to the ridge, with the leeward 
eastern side consisting of a steep and collapsed the 
topography scraped by avalanches (Fig. 9A), with narrow 
ridges with sheer sides sporadically present. On the very 
steep and collapsed slope, no grasses grow and even no 
transitional grasslands of Poa or Carex grass vegetation 
can be found. Only in the area around the summit of Mt. 
Rengedake there exists a topography of wide and gentle 
gravel areas with large colonies of Dicentra peregrina 
with very few Pinus pumila (Fig. 9F). These factors suggest 
that the distribution of alpine butterflies, O. norna, E. 
ligea, and E. neriene that feed on Poa and/or Carex grasses, 
is also restricted by the geography and vegetation around 
the Harinoki Gap. Another explanation of the vicarious 
distribution of O. norna in the Harinoki Gap has been 
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proposed, in which it was caused by the sedimentation of 
volcanic products of Mt. Tateyama approximately 100,000 
years ago (Itoh, 2010). However, the distribution patterns 
es oI of these three alpine butterflies strongly suggest that 
SES 6 % 2 1 3 species-specific ecological factors were the major cause. 
on) | | | 

= Cie aa The process of range expansion and contraction of alpine 
SO FÈ organisms due to climatic changes in the Quaternary period 

has been well known as a north-south migration in the 
Japanese archipelago (Ono, 1990). However, it could be 
possible that more than one refugium existed in the Japanese 
archipelago because of its complex mountain topography 
and the resulting micrometeorological environments. As 
a result of isolation by taking refuge in more than one 
refugium during the warm periods, genetically differentiated 
populations have survived till now as multiple lineages 
by expanding and contracting their distribution areas in 
adaptation to later glacial cycles (the multiple refugia 
model). The first example of such a multiple refugia model 
was suggested in the study on an alpine butterfly Erebia 
neriene (Nakatani et al., 2007a). The present study also 
suggests the multiple refugia model for another alpine 
butterfly E. ligea, and it is further suggested that more 
organisms may well be subject to the same considerations. 
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Fig. 8. Vertical distributions of three alpine butterflies in the Harinoki gap, in the central Hida Mts. 
Solid circle: Erebia neriene, open circle: Erebia ligea, star: Oeneis norna. 
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Fig. 9. Landscapes around the Harinoki gap, in the central Hida Mts. 


A: Steep ridge near Mt. Fudodake. East (left) side with collapsed cliff, west side covered with tall forest, 
with no grassland. 

B: Rocky cliff, between Mt. Subaridake and Mt. Harinokidake. Sporadic dried herbaceous communities, 
no Erebia species. 

C: Mt. Akazawadake overlooked from Mt. Subaridake. Some alpine grassland spreads on the leeward 
slope (right, east side). Erebia species and Colias palaeno inhabit. 

D: Mt. Subaridake (left) to Mt. Akazawadake and Mt. Narusawadake (right, NNE). Light green zones 
near the ridge are grasslands where Erebia species are found. 

E: Saddle point area, Nanakura-nokkoshi (right, east side). Erebia species found. 

F: Wide and gentle ridge near the summit of Mt. Rengedake. Dicentra peregrina (pink flowers) abundant 
and some Carex. No alpine butterflies. 
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摘 要 


第 四 紀 日 本 列島 に お ける クモ マ ベ ニ ヒ カ ゲ の 
分 布 変遷 史 (中 谷 貴 毒 ・ 宇 佐美 真一 ・ 伊 藤 建 夫 ) 


日 本 列島 は 氷河 期 に は 大 陸 と 繋が っ て いた 時 代 が あり , そ 
の 期間 に は 大 陸 起 源 の 北方 系 生物 は 分 布 を 拡大 し , 間 氷 期 
に は これ ら の 生物 は 大 陸 へ 避難 し , ある い は 一 部 の 個体 群 
を 除い て 絶滅 し た と 考え られ て きた . LL, 筆者 ら の 高 
山 虹 に 関す る 研究 に よる と , ベニ ヒ カ ゲ は 古い 時 代 の 水 河 
期 に 渡来 し , 間 氷 期 に も 日 本 列島 の 高山 常 で 複数 の 個体 群 
が 生き 残り , 互い に 生殖 隔離 され た 結果 列島 内 で 複数 の 系 
統 に 分 化し た 後に , 次 の 氷河期 に 分 布 を 拡大 する と いう サ 
イク ル を 繰り 返し て きた こと が 明らか に な っ た (複数 レフ 
ュ ジ ア ・ モ デル ). これ は ヨ ツ バシ オガ マ 他 数 種 の 高山 植 
物 で 明らか に され た シナ リオ , 初め の 氷河期 に 大 陸 か ら 渡 
来 し た 系 統 が 間 氷 期 に 本 州 中 部 山岳 で 生き 残り , 次 の 氷河 
期 に 新た に 侵入 し た 系 統 が 東北 地方 以 北 に 分 布 す る と いう 
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時 間 差 侵入 仮説 ( 単 二 レフ ュ ジ ア ・ モ デル ) と は 異な る . 


本 研究 で は サハ リン を 含む 日 本 列島 か ら 17 の ハ プ ロ タ イ 
プ を 見 出し た (Table 1). 複数 地域 また は 複数 サン プル か 
ら 検 出さ れ た 系 統 的 意義 を 有する ハプロタイプ は , 本 州 で 
は 飛騨 山脈 北部 ・ 白 山 , 飛騨 山脈 南部 , 八ヶ岳 , 木曽 山脈 , 
赤石 山脈 の 5 系統 に , また 北海 道 で は 大 雪山 の 高 標高 部 と 
山 薦 の 低 標高 部 の 2 系 統 に 分 か れ て いる (利尻 島 高 山 帯 に 
も 孤立 し た 個体 群 が 分 布 す る が 未 検 , Fig. 3). サハ リン . dt 
海道 , 本 州 の 個体 群 は , 過去 の 異 所 的 分 断 に より 分 断 分 布 
を 成 し て いる 事 が 明らか で ある . 本 州 で は 飛騨 山脈 で 複数 
の ハプロタイプ が 混 生 し て お り (Fig.4), 複数 の イベ ント 
に よっ て 現在 の 分 布 が 形成 され た と 考え られ る の で , NCPA 
に より 過去 の 分 布 変遷 史 を 推定 し た 結果 , 5 つの クレ ー ド 
で 統計 的 に 有意 な イベ ント が 推定 され た (Table 2). それ 
に よる と , クモ マ ベ ニ ヒ カ ゲ の 中 部 山岳 地域 に お ける 分 布 
の 変遷 は . 分 断 と 拡散 の 繰り 返し で ある こと が 示唆 され た . 


日 本 列島 に お ける 分 布 変遷 


日 本 列島 へ 進出 し た クモ マ ベ ニ ヒ カ ゲ の 個体 群 は . サハ リ 
ン ・ 北 海道 ・ 本 州 の 現在 の 分 布 域 を 包含 する 地域 に 分 布 を 
広げ た . その 後 の 温 暖 期 に サハ リン , 北海 道 , 本 州 の レフ 
ュ ジ ア に 分 断 さ れ , それ ぞ れ が 別々 の 系 統 に 分 化し た . 次 
の 温暖 期 に , 本 州 の 中 部 山岳 地域 で は 飛騨 山脈 北部 系 統 , 
同 南部 系 統 , 赤石 山脈 ・ 木 曽山 脈 系 統 の 3 系 統 に 分 断 さ れ 
た . 続く 氷河 期 に 分 布 を 拡大 し , 飛騨 山脈 の 北部 系 統 と 南 
部 系 統 は , 後 立山 連峰 の 針 ノ 木 岳 付近 (以後 針 ノ 木 ギ ャ ッ 
プ と 呼ぶ ) で 混 生 地帯 を 形成 し た . さら に その 後 の 温暖 期 
に 現在 見 る よう な 離散 分 布 が 形成 され た . 


ベニ ヒ カ ゲ と の 比較 に よる 系 統 地 理 的 な 特徴 


広域 に 分 布 す る ハプロタイプ の 系 統 関係 を 概観 する と (Fig. 
7), 両 種 と も に 飛騨 山脈 と 赤石 山脈 産 の ハプロタイプ の 間 
に 大 き な 遺 伝 的 差異 の ある こと が 示さ れる . 初期 の 単 一 な 
遺伝 的 組成 を も つ 集 団 が , その 後 の 温暖 期 に 分 布 を 縮小 す 
る 過程 で 飛騨 山脈 と 赤石 山脈 に 分 布 す る 二 つ の 個体 群 に 
分 断 さ れ た 結果 , 二 つ の 系 統 に 分 岐 し た こと を 示唆 し て い 
る . 飛騨 山脈 と 赤石 山脈 の レフ ュ ジ ア に 源 を 発する 2 系 統 
の 遺伝 的 距離 は , ベニ ヒ カ ゲ と クモ マ ベ ニ ヒ カ ゲ と の 間 で 
差異 が あり , 分 岐 年 代 に は 若干 の 差 が あっ た と 老 え られ る . 
また ベニ ヒ カ ゲ と クモ マ ベ ニ ヒ カ ゲ 共 に ., 赤石 山脈 の 系 統 
は , 飛騨 山脈 系 統 の 北部 集団 と より 近 縁 で ある 事 は , その 
後 の 気 候 変動 に 適応 し て 分 布 を 拡大 し た ルー ト が 両 種 で 類 
似 し て いる 可能 性 を 示唆 し て いる . 


飛騨 山脈 で は , ベニ ヒ カ ゲ は 単 一 の ハプロタイプ が 産 す る 
が , クモ マ ベ ニ ヒ カ ゲ で は 針 ノ 木 ギ ャ ッ プ で 混 生 地 を 挟ん 
で 南北 2 系 統 に 分 か れ て お り , 両 系 統 の 間 に 2 塩基 の 差 が 
認め られ る 点 が 大 きく 異な る . 白山 山系 は 飛騨 山脈 と は 地 
理 的 に 非常 に 離れ て いる が ベニ ヒ カ ゲ , クモ マ ベ ニ ヒ カ ゲ 
(クモ マ ベ ニ は 飛騨 山脈 の 北部 系 統 ) 共に 同じ ハ プ ロ タ イ 
プ が 分 布 し て お り , 二 つ の 山系 の 個体 群 間 で は きわ め て 最 
近 ま で 遺伝 的 交流 の あっ た こと が 示唆 され る . 八ヶ岳 で は , 
ベニ ヒ カ ゲ は 飛騨 山脈 と 同一 の ハプロタイプ が , また クモ 


マ ベ ニ ヒ カ ゲ は 飛騨 山脈 の 南部 系 統 と 1 塩基 差 の 近 縁 な ハ 
プロ タイ プ が 見 出さ れ , 両 地域 の 個体 群 の 近 縁 性 が 示唆 さ 
れる . これ に 対し て 木曽 山脈 で は . ベニ ヒ カ ゲ は 飛騨 山脈 
と 共通 の ハプロタイプ が 見 出さ れる の に 対し て , クモ マ ベ 
ニ ヒ カ ゲ で は 赤石 山脈 と 近 縁 で あり , 木骨 山脈 の 個体 群 の 
分 布 変遷 は 両 種 の 間 で 異な っ て いた 可能 性 が 示唆 され る . 
一 方 サハ リン と 北海 道 の 個 体 群 に 関し て は , ベニ ヒ カ ゲ は 
共通 の ハプロタイプ が 分 布 す る な ど , 最近 まで 遺伝 的 交流 
が あっ た こと が 示唆 きれ る が , クモ マ ベ ニ ヒ カ ゲ の 個体 群 
は 遺伝 的 に 非常 に 異な っ て お り , 古く か ら 生 殆 隔 離 が 続い 
て いる こと が 示唆 され た . 


ベニ ヒ カ ゲ と クモ マ ベ ニ ヒ カ ゲ は 同じ 属 に 含ま れる 近 緑 種 
で あり , また 生息 環境 も 似 て いる が , いく つか の 地域 で は 
異な る 分 布 変遷 更 を た どっ た よう だ . この よう に 種 に よっ 
て 生殖 隔離 の 始ま っ た 時 期 や 場所 が 異な る 事例 は ヨー ロッ 
パ で も 知ら れ て いる . Erebia medusa は ヨー ロッ パ で は 針 葉 
樹林 内 の 草原 を 主たる 生息 地 と し て いる が . ルー マニ ア と 
ブル ガリ ア の 28 集 団 に つい て 調べ た 研究 に よる と . ドナ 
ウ 河 を 挟ん で 南北 2 系 統 に 分 断 さ れ て お り , 最終 氷 期 に そ 
れ ぞ れ の 集団 が 分 布 を 拡大 し た も の の ドナ ウ 河 を 越え て 遺 
伝 的 交流 が 成 き され る こと は な か っ た と し て いる . 一 方 E. 
medusa と 同様 に 針葉樹 林内 の 草原 を 主たる 生息 地 と し て 
広域 分 布 す る Erebia euryale で は , ルー マニ ア 産 と ブル ガリ 
ア 産 は 遺伝 的 に よく 似 て お り , 最終 水 期 以降 に ドナ ウ 河 を 
越え て 遺伝 的 交流 が あっ た と する 研究 が ある . この よう に 
大 陸 内 で し か も 生息 環境 が 類似 し た 種 で も , 第 四 紀 の 水 河 
サイ クル に 対す る 適応 は 種 に よっ て 異な る 事例 が ある よう 
に , 日 本 列島 内 に お ける 第 四 紀 の 気候 変動 に 対す る 高山 蝶 
の 適応 は . 種 特異 的 な いろ いろ な パタ ー ン が 存在 する こと 
が 強く 示唆 され る . 


針 ノ 木 ギャ ッ プ の 生態 的 意義 


すでに 述べ た よう に , 飛騨 山脈 に お ける クモ マ ベ ニ ヒ カ ゲ 
の ハプロタイプ は , 鹿島 槍ヶ岳 か ら 鳥 帽子 岳 に 至る 針 ノ 木 
ギャ ッ プ で 2 系 統 の 混 生 地 淀 が 見 られ る (Fig.4). タカ ネ 
ヒ カ ゲ は 標高 約 2.700 m AEO RETI DDA TYRD 
け て の 高山 帯 に の み 生 息 す る 真性 高山 将 で , 針 ノ 木 ギ ャ ッ 
プ 付 近 で は 爺 ヶ 后 か ら 鳥 帽子 項 の 間 で 分 布 を 欠い て お り , 
雪 倉 后 か ら 鹿 島 槍 ヶ 后 ・ 布 引山 に か け て 分 布 する 北部 系 統 
と , 烏帽子 岳 以 南 に 分 布 する 南部 系 統 に 分 岐 し て いる . 一 
方 , ベニ ヒ カ ゲ は 針 ノ 木 ギ ャ ッ プ を 含む 飛騨 山脈 全域 に 単 
一 の 系 統 が 分 布 す る . 針 ノ 木 ギ ャ ッ プ 付近 に お ける これ ら 
3 種 の 高山 蝶 の 分 布 状況 は , 種 の 標高 に 対す る 適応 の 度合 
い を 反映 し た も の と な っ て いる (Fig.8). 針 ノ 木 ギ ャ ッ プ 
付近 の 地形 お よび 植生 を みる と (Fig.9), 全体 に 標高 が 約 
2500-2600 m č Rozo CEROS l DBM CH DIT 
BY), 冬季 季節 風 の 風 上 に 当たる 尾根 の 西側 で は 樹林 が 尾 
根 ま で 迫り , 風下 側 の 東側 は 雪 月 に よっ て 削ら れ た 急 唆 な 
崩落 地形 と な り , ある い は 両側 の 切り 立っ た 狭い 尾根 が 断 
続 的 に 見 られ る . AWE HAMIL ABP AVY TH 
科 の 乱 移 途上 の 草地 すら 見 られ な い . 蓮華 后 の 頂上 付近 に 
の み 広 い 緩 傾斜 の 砂 秦 帯 が あり , コマ クサ が 大 群落 を 形成 
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する が ハイ マツ は あま り 生 えて いな い . この よう な 植生 が , 
イネ 科 や カヤ ツ ッ リ グ サ 科 植 物 を 食 草 と する 高山 上 嘩 た ち ( タ 
カネ ヒ カ ゲ , クモ マ ベ ニ ヒ カ ゲ , ベニ ヒ カ ゲ ) の 分 布 を 規 
fl し て いる も の と 考え られ る . 人 針 ノ 木 ギ キャップ 付近 に お け 
る タカ ネ ヒ カゲ の 不 連 続 分 布 の 原因 を 約 10 万 年 前 の 立山 
噴火 に よる 噴出 物 の 堆積 に 求め る 説 も ある が , 3 種 の 高山 
虹 に み ら れ る 分 布 パ ター ン は , それ ぞ れ の 種 が 持つ 高度 適 
応力 の 強 さ を 反映 し て お り , 生態 的 要因 が より 強く 働い た 
結果 で ある と 考え られ る . 


従来 は 日 本 列島 の 生物 相 形 成 に 関し て , 大 陸 か ら の 複数 回 


の 進出 に よっ て 氷河 サイ クル に 適応 し た と する 事例 ( 単 一 
レフ ュ ジ ア ・ モ デル ) が 指摘 きれ る ケー ス が 多かっ た . し 
か し 筆者 ら に よる ベニ ヒ カ ゲ や タカ ネ ヒ カゲ , さら に は 今 
回 報告 し た クモ マ ベ ニ ヒ カ ゲ の 研究 で 示唆 きれ る よう に , 
高山 性 生物 の 種 に よっ て は 古く か ら 日 本 列島 に 侵 出し , 日 
本 列島 内 で 氷河 サイ クル を 通じ て 分 布 の 分 断 ・ 拡 張 を 繰 返 
し て きた ケー ス (複数 レフ ュ ジ ア ・ モ デル ) が 少な く な い 
こと が 伺わ れる . 


(Received April 23, 2012. Accepted November 8, 2012) 
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